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The static conductivity of the dilute magnetic semiconductor Gai-xMn^As is calculated using the 
memory function formalism and time-dependent density-functional theory to account for impurity 
scattering and to treat Hartree and exchange interactions within the hole gas. We find that the 
Coulomb scattering off the charged impurities alone is not sufficient to explain the experimentally 
observed drop in resistivity below the ferromagnetic transition temperature: the often overlooked 
scattering off the fluctuations of localized spins is shown to play a significant role. 
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The perspective of utilizing the charge and the spin 
of the electrons for new electronic device applicati- 
ons has generated tremendous interest in the field of 
spintronics.^ A unique combination of magnetic and 
semiconducting properties makes dilute magnetic semi- 
conductors (DMSs) very attractive for various spintronics 
applications.^ Among the family of DMSs, much atten- 
tion has been paid to Gai-ajMn^; As since the discovery of 
its relatively high ferromagnetic transition temperature,^ 
with a current record of Tc = 185 K.^ 

Unlike most other III-V DMSs, the nature of the itin- 
erant carriers in Gai-a^Mn^^As is still under debate. It 
is widely accepted that for low-doped insulating samples 
the Fermi energy lies in a narrow impurity band. For 
more heavily doped, high-Tc metallic samples there are 
strong indications that the impurity band merges with 
the host semiconductor valence band forming mostly 
host-like states at the Fermi energy with some low-energy 
tail of disorder-related localized states.^ First-principles 
calculations^i^ii have so far not been fully conclusive re- 
garding the nature of the itinerant carriers in this case, 
and further theoretical studies continue to be necessary. 
The question is thus, in essence, whether the valence 
band^ or impurity band^ picture is more adequate to de- 
scribe the various experimental results in Gai-ajMn^^As. 

The purpose of this paper is to present a study 
which supports the valence band picture for electronic 
transport properties and for the optical conductivity in 
Gai-a^Mn^^As. In this material, unlike in II- VI DMSs, 
the magnetic ions in substitutional positions act as ac- 
ceptors delivering holes and producing not only localized 
spins but also charged defects. We will argue that it is im- 
portant to treat disorder and many-body effects beyond 
the simple relaxation time and static screening models 
which were used in earlier theoretical studies. ^^^^^^^^^^^ 

Lopez-Sancho and Bre}^ considered the temperature 
dependence of the Coulomb scattering off the accep- 
tor centers and found that the carrier relaxation rate 
is reduced by around 20% in the ferromagnetic phase, 
consistent with the experimentally observed drop in 
resistivity^i^ii^ii^ This drop was attributed entirely to 
the effects associated with the scattering off Coulomb 
disorder. The main mechanism found to be responsi- 



ble for the observed drop in resistivity was the change 
of the semiconductor band structure: during the transi- 
tion from the paramagnetic to the ferromagnetic state, 
the giant spin splitting of the energy bands substantially 
modifies the shape of the Fermi surface, thus altering the 
possible scattering wave vectors and, consequently, the 
magnitude of scattering matrix elements. 

These findings speak in favor of the valence band pic- 
ture of Gai-a^Mn^^As. However, the model of Ref. I12I 
employed a simplified treatment for the screening of the 
charge disorder by itinerant carriers, neglecting the ex- 
change part of the electron-electron interaction within 
the hole liquid. Furthermore, the scattering off the fiuc- 
tuations of localized spins was ignored. We will show that 
both effects play an important role in spin polarized sys- 
tems and should be included in the valence band picture 
model of the itinerant holes in Gai-a^Mn^^As. In fact, 
our calculations suggest that the previously suggested 
origin^^ of the resistivity drop in the ferromagnetic phase 
should be revised: the main reduction of the scattering 
rate comes from the suppression of the fiuctuations of 
localized spins in the magnetically ordered state. 

Earlier we developed a theory of transport in charge 
and spin disordered media which combines a multiband 
k • p approach for an accurate description of the valence 
band states with a more comprehensive treatment of dis- 
order and electron-electron interaction. ^^^^^ Our theory is 
based on an equation of motion approach for the current- 
current response functio n^^i^^ and has some similarity 
with the memory function formalism 

In this paper we apply our formalism to describe the 
transport properties of Gai-a^Mn^^As in the static regime. 
Specifically, we focus on the pronounced drop in resis- 
tivity below Tc which has been observed for optimally 
annealed metallic samples. To investigate this phe- 
nomenon let us first look at the standard expression for 
the static conductivity obtained from the semiclassical 
Boltzmann equation: 
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where the summation is over the wave vector k and the 
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energy band index n. The part of Eq. ([T]) which is most 
sensitive to temperature is the carrier scattering rate t~^. 
The task is therefore to derive a microscopic expression 
for r~J^ which accounts for ah relevant scattering mech- 
anisms, as weh as for electronic many-body effects. 

Our model for the itinerant carriers in Gai-a^Mua^As is 
that of a system with charge and spin disorder described 
by the Hamiltonian 



limit, which suffices to illustrate the general structure: 
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where the four-component disorder scattering potential 



-T 



(3) 



is coupled to the four-component charge and spin density 
operator of the band carriers. 
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Here, (/i = 1, +, — , is defined via the Pauli matrices, 
where is the 2x2 unit matrix, = (cr^ ± i<j^)/2, 
and |iin,q) are the two-component Bloch function spinors 
with wave vector q and band index n. The summation 
in Eq. (|3]) is performed over all defects. 

The general case of multiple types of defects, includ- 
ing defect correlations, was considered in Ref. 17. For 
simplicity we here include only the most important de- 
fect type, namely randomly distributed manganese ions 
in gallium substitutional positions (Mnca)- Our model 
treats localized spins as quantum mechanical operators 
coupled to the band carriers via a contact Heisenberg in- 
teraction featuring a momentum-independent exchange 
constant J. We use the value of VJ = 55 meVnm^, 
which corresponds to the widely used DMS ^—d exchange 
constant A^o/^ = 1.2 eV.^ The z-axis is chosen along the 
direction of the macroscopic magnetization. 

We obtain the following expression for the tensor of 
Drude-like frequency and momentum dependent relax- 
ation rates: 

= ^Z^E^"^/5(^A.(-k)W.(k)) (5) 
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where a^jS = x^y^z are Cartesian coordinates, n is the 
carrier concentration, and Xp/^p^(k, c^) are charge- and 
spin-density response functions associated with the oper- 
ators dH). The superscript c in Eq. ([5]) refers to a clean 
(defect-free) system, and r^^ stands for those additional 
contributions which arise in magnetically ordered sys- 
tems only. Here, we just give the clean-system and q = 
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The complete expression for r^^ as well as details of the 
derivation of Eq. ([5]) can be found in Ref. [H. 

All information about the itinerant carriers, includ- 
ing band structure and electron-electron interaction, is 
contained within the set of charge- and spin-density re- 
sponse functions Xp/^p^ (k, oo). Strictly speaking, these re- 
sponse functions correspond to the disordered system, 
and Eq. ([5]) should be calculated self-consistently.^^ We 
here assume that the disorder is sufficiently weak so that 
we can approximate ([5]) by expanding to second order in 
the disorder potential Z^(k), and thus replace Xp/^p^ (k, oo) 
by its clean system counterpart Xp^p^- (k, a;). 

To account for the complexity of the band structure, 
we use a standard 8-band k • p approach with contribu- 
tions from the remote bands taken up to the second order 
in the wave vector The mean-field part of the p — d ex- 
change interaction between itinerant holes and localized 
spins causes a spin splitting of the bands of the semicon- 
ductor host material. Technical details of this multiband 
linear-response approach will be published elsewhere. 

The major advantage of Eq. ([5]) is that it can be com- 
bined in a straightforward manner with time-dependent 
density functional theory (TDDFT),^^ which allows us 
to describe electron-electron interaction effects, including 
correlations and collective modes, in principle exactly. In 
TDDFT the charge- and spin-density response functions 
of the interacting system are written as follows:^^ 
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where xo denotes the matrix of response functions of the 

noninteracting system, ^(g) is the Hartree part of the 

electron-electron interactions, and fxc is the matrix of 

exchange and correlation kernels. All quantities in Eq. 
([7j) are 4x4 matrices; according to Eq. (|3]), the first 
component is charge, and the other components are spin 
+, — , and As a simplification we use only the exchange 
part of fxc and apply the adiabatic local spin density 

approximation. The local field factors for partially spin 
polarized system were calculated according to Ref. |28|. 

Since the mean-field part of the p — d exchange interac- 
tion is extracted from the disorder Hamiltonian, the total 
relaxation rate ([5]) can be separated into contributions 
associated with Coulomb disorder and with fiuctuations 
of the localized spins. The transverse component (per- 
pendicular to the magnetization) of the relaxation rate 
tensor in the long wave length (q = 0) and static {uo 0) 
limit then has the form 
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FIG. 1: Temperature dependence of the carrier relaxation 
rates in Gao.95Mno.05As. The vertical line indicates Tc. (a) re- 
laxation rates associated with the scattering off Coulomb dis- 
order, calculated within different screening models (see text 
for details), (b) total relaxation rate (solid line) and contribu- 
tions from scattering off Coulomb disorder (dotted line) and 
off localized spin fluctuations (dashed line) . 



where the charge disorder contribution is 
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and the contribution from the fluctuations of the local- 
ized spins is given by 
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Here, rii denotes the concentration of Mnca defects. U {k) 
represents the Coulomb potential of a single acceptor 
center screened by the dielectric constant of the host 
material, where we take e = 13 for GaAs; the screen- 
ing by the electron liquid is absorbed in the band-carrier 
response functions. The angular brackets in Eq. (p!Q|) 
denote the thermodynamic average with respect to the 
magnetic subsystem Hamiltonian Hm in Eq. ([5]). We as- 
sume Hm to be a pairwise Heisenberg-like Hamiltonian. 
In our calculations we use the experimental value of Tc 
as an input parameter and apply the standard mean field 
approach to obtain the temperature dependence of ther- 
modynamically averaged quantities in Eq. (p!Q|) . 

Fig. [1] presents the typical temperature dependence 
of the carrier relaxation rate in Gai-a^Mn^^As obtained 
within our model. The calculations were performed for 



a system with Mn concentration x = 0.05 and carrier 
concentration of p = 0.3 holes per Mn in Ga substitu- 
tional positions. The left panel shows the contributions 
from the scattering off Coulomb disorder with screen- 
ing effects accounted for by three different methods. The 
dotted line corresponds to screening described within the 
Thomas- Fermi approximation for paramagnetic systems. 
The screening here is temperature independent and the 
20% drop in the scattering rate in the ferromagnetic 
phase is entirely due to modification of the possible scat- 
tering wave vectors, the mechanism described in Ref ITU. 

If, however, we allow the change of the band structure 
to affect the screening as well, e.g. on the random-phase 
approximation (RPA) level (first two terms in Eq. ([7j) 
and dashed line in Fig. [T^), the drop in the resistivity 
in the ferromagnetic phase is significantly reduced. This 
effect was also considered in Ref. UH. But if we now go 
further and include the exchange part of electron-electron 
interaction in Eq. ([7j), then the drop in the resistivity 
is completely washed out, see the solid line in Fig. [T^. 
Moreover, for some parameters, the trend is reversed and 
the scattering off Coulomb disorder actually increases in 
the ferromagnetic phase. 

The explanation for this behavior lies in the exchange 
part, which counteracts the larger Hartree part of the 
electron-electron interaction. It reduces, therefore, the 
screening of the Coulomb disorder potential calculated 
within RPA. Therefore we have an overall increase in the 
charge relaxation rate once the exchange part of electron- 
electron interaction is taken into account, see Fig. [1^. 
On the other hand, the exchange part of the electron- 
electron interaction is more pronounced for spin polarized 
systems, resulting in a stronger reduction of the screening 
of the disorder potential and thus causing an increase 
of the Coulomb scattering in the ferromagnetic phase. 
This process, neglected in Ref.fHI, competes with and, for 
some parameters, reverses the reduction of the relaxation 
rate due to band-structure related modifications of the 
scattering wave vectors. 

It is thus apparent that the scattering off the Coulomb 
disorder potential alone cannot be responsible for the ex- 
perimentally observed drop in resistivity. The other pos- 
sible contribution is the scattering off the fluctuations 
of localized spins. In Fig. [TJd we plot the temperature 
dependence of the scattering rate for both mechanisms. 
The scattering off spin fluctuations is often overlooked 
since its magnitude is substantially smaller than that of 
the Coulomb scattering. Due to effective suppression of 
spin fluctuations in the ferromagnetic phase, however, the 
temperature dependence of this relaxation mechanism is 
much more pronounced. Indeed, in a fully spin polarized 
state, the scattering takes place only off the quantum 
fluctuations of localized spins. The total relaxation rate 
dHI), which is the sum of both contributions, restores its 
20% drop during transition from paramagnetic to ferro- 
magnetic phase. The majority of this drop is found to 
be due to the suppression of the scattering off localized 
spin fluctuations. 
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FIG. 2: Static resistivity of Gao.95Mno.05As normalized to the 
paramagnetic state value. Solid and dashed lines: compensa- 
tion levels p = 0.3 and p = 0.5 (number of holes per Mnca). 
The vertical line indicates Tc — 110 K. Experimental points 
are from Ref. .14 . 

In Fig. [2] we plot the temperature dependence of the 
static resistivity of Gao.95Mno.05As, normalized to the 
paramagnetic state value. The calculation was done ac- 
cording to Eqs. ([1]) and (|8|)- (p!Q|) . Solid and dashed lines 
correspond to different levels of compensation in the sys- 
tem, 0.3 and 0.5 hole per substitutional Mn, respectively 
(in practice, this number is difficult to control). The open 
squares represent the experimental data of Ref. [13. The 



theory demonstrates good agreement with experiment. 

In summary, we have developed a theory of transport 
in spin and charge disordered media within the valence 
band picture of metallic GaMnAs. The approach com- 
bines the multiband k • p description of the semiconduc- 
tor band structure with a microscopic treatment of dis- 
order and dynamical electron-electron interaction by the 
methods of the memory-function formalism and TDDFT. 
We applied our formalism to describe the experimentally 
observed drop in static resistivity of GaMnAs in the fer- 
romagnetic phase. 

This problem had been addressed before in Ref. 12, but 
with a model that was lacking some important features 
such as scattering off the fluctuations of localized spins 
and electron-electron interactions beyond RPA. Similar 
to Ref. 12, we obtained agreement with the experimental 
observations, but the underlying physics is quite differ- 
ent. Much of the drop of resistivity in the ferromagnetic 
phase is found to be due to the suppression of localized 
spin fluctuations in the magnetically ordered state. 

To conclude, we have developed a theoretical descrip- 
tion of itinerant carriers in DMSs within the valence band 
picture that accounts for band structure, scattering from 
Coulomb and magnetic impurities, and screening via dy- 
namical many-body effects. An accurate description of 
static transport properties in GaMnAs involves a subtle 
interplay of all these ingredients. Our approach is also 
suitable for frequency-dependent properties such as the 
optical conductivity, as will be discussed elsewhere. 

This work was supported by DOE Grant No. DE- 
FG02-05ER46213. 
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